Submonolayer epitaxy with impurities 



Miroslav Kotrla 

Institute of Physics, Czech Academy of Sciences, Na Slovance 2, 182 21 Praha 8, Czech Republic 

Joachim Krug 

Fachbereich Physik, Umversitdt GH Essen, 45117 Essen, Germany 
Pavel Smilauer 

Institute of Physics, Czech Academy of Sciences, Gukrovarnickd 10, 162 53 Praha 6, Czech Republic 
(Submitted to Physical Review B on January 25, 2000) 

The effect of impurities on epitaxial growth in the submonolayer regime is studied using kinetic 
Monte Carlo simulations of a two-species solid-on-solid growth model. Both species are mobile, 
and attractive interactions among adatoms and between adatoms and impurities are incorporated. 
Impurities can be codeposited with the growing material or predeposited prior to growth. The 
activated exchange of impurities and adatoms is identified as the key kinetic process in the formation 
of a growth morphology in which the impurities decorate the island edges. The dependence of the 
island density on flux and coverage is studied in detail. The impurities strongly increase the island 
density without appreciably changing its power-law dependence on flux, apart from a saturation of 
the flux dependence at high fluxes and low coverages. A simple analytic theory taking into account 
only the dependence of the adatom diffusion constant on impurity coverage is shown to provide 
semi-quantitative agreement with many features observed in the simulations. 



I. INTRODUCTION 

Recent progress in the fabrication of atomically smooth 
interfaces by Molecular Beam Epitaxy (MBE) has lead 
to an increasing appreciation of the dramatic, detrimen- 
tal or beneficial, effects that small amounts of impu- 
rities may have on the morphology of growing films. 
Adsorbates acting as surfactants can stabilize layer-by- 
layer growth of metal"'^"^ and semiconductor'"'^ surfaces. 
On the other hand, for the simple case of Pt(lll) ho- 
moepitaxy it was recently shown^ that minute cover- 
ages of CO strongly increase the step edge barriers for 
interlayer transport, thus enhancing three-dimensional 
mound growth.^" The effect of additional surface species 
on growth and nucleation is of obvious importance also in 
more complex, technologically relevant deposition tech- 
niques such as chemical vapor deposition. In either case 
the detailed atomistic kinetics and energetics of the in- 
teraction between adsorbate and deposited material in- 
fluence the growth mode to a degree which makes it very 
difficult to formulate general rules for large classes of 
growth systems. 

As a first step towards an improved understanding of 
the generic effects of impurities on epitaxial growth, in 
the present paper we introduce a minimal model which, 
we hope, is simple enough to extract some insights of 
fairly general validity, and yet possesses sufficient fiex- 
ibility to include most physically relevant microscopic 
features. The model is based on the standard solid- 
on-solid description for the growth of a simple cubic 
crystal. ^^'"'^'^ The impurities are represented by a second 
particle species, which can be codeposited with the grow- 
ing material or predeposited prior to growth. Impurities 



diffuse and interact attractively with the deposit atoms 
(adatoms) , but they do not attract each other and hence 
do not nucleate islands. The details of the model are de- 
scribed in Section II. A brief account of some preliminary 
results was given in an earlier communication.^** 

Despite its simplicity, the model contains a large num- 
ber of parameters: impurity and growth fiuxes, substrate 
temperature, and energy barriers for half a dozen ki- 
netic processes. To focus our efforts, we concentrate 
on modeling a situation in which the impurities decorate 
the island edges, forming a monatomic chain along the 
island perimeter. Preferential adsorption of impurities 
at step edges is suggested by bond counting arguments, 
and has often been invoked to explain the strong effect 
of submonolayer adsorbate coverages on growth behav- 
ior, e.g., through a change of the barrier for interlayer 
transport. ^'^'^^ 

It will be shown in Section III that the growth of dec- 
orated islands requires, in addition to a suitable choice 
of binding energies, the possibility of impurity-adatom 
exchange. Such a process, which is a two-dimensional 
analog of the exchange mechanism responsible for the 
fioating of surfactants in multilayer growth,^ was recently 
demonstrated to play a crucial role in the submonolayer 
homoepitaxy of Si(OOl) in the presence of hydrogen. -'^^ 

In Section IV the influence of the adsorbates on the 
island density is investigated, both for codeposited and 
predeposited impurities. In the absence of impurities the 
scaling relation 

N - {F/DY (1) 

between island density N , deposition flux F and adatom 
diffusion coefficient 
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D is well established theoretically,^^ numerically^''^^'' 
and experimentally.^'* In two dimensions, rate equation 
analysis^^"^^ yields the expression 



for the exponent x iii terms of the size i* of the largest 
unstable cluster. 

There are several conceivable mechanisms by which 
attractive impurities^^ could alter the relationship (1). 
First, impurities may act as nucleation centers, thus ef- 
fectively decreasing i* and therefore in the extreme 
case of immobile adatom traps the limit of spontaneous 
nucleation with i* = x = would be realized. ^^'^^ Sec- 
ond, impurities decorating the island edges may induce 
energy barriers to attachment. Kandel^'^ has shown that, 
provided these barriers are sufficiently strong, the expo- 
nent X in (1) is increased such that (2) is replaced by 
X = 2z*/ {i* + 3). Both mechanisms imply an increase of 
the island density compared to the case of pure homoepi- 
taxy. 

Our simulations indicate that none of these two mech- 
anisms are operative under the conditions used in our 
model: The addition of impurities is found to increase 
the island density in all cases, but the scaling of N 
with the flux F remains unaffected within the accuracy 
of the simulation. An analysis of the relevant micro- 
scopic processes^* reveals that, within our model, even 
completely decorated island edges do not provide effi- 
cient barriers to attachment, and therefore the scenario 
of Kandcl^^ does not apply. 

In view of (1), an increase of the island density at fixed 
X suggests that the main effect of the impurities is to re- 
duce the mobility D of adatoms. The reduction of the 
adatom mobility has been identified as the most impor- 
tant mechanism contributing to the surfactant action of 
Sb on Ag.^'^'^* Decreasing D reduces the island size and 
favors the growth of ramified, rather than compact is- 
lands. Both effects enhance intcrlayer transport, since 
the adatoms landing in the second layer have more op- 
portunities to descend, and thus promote layer-by-layer 
growth. In our model the island size decreases but the 
island shapes remain compact, because the edge decora- 
tion facilitates edge diffusion (see Section III). 

For an analytic description of the relation between im- 
purity coverage, adatom mobility and island density, in 
Section V we develop a simple rate equation approach 
which provides a semi-quantitative explanation for many 
(though not all) features observed in the simulations. 
Some conclusions and open questions are formulated in 
Section VI. 



II. MODEL 

The growth model employed in this work has been 
briefly described in an earlier paper. It is a solid-on- 
solid model with two surface species A and B, where 



A-particlcs correspond to the growing material, and B- 
particlcs represent the impurities. The simulation starts 
on a flat substrate composed only of y4-atoms. The; ba- 
sic microscopic processes are deposition and migration; 
desorption is not allowed. Two deposition modes are 
considered: (i) simultaneous deposition (codeposition) of 
both species and (ii) predeposition of a certain impurity 
coverage prior to growth. In the case of codeposition the 
fluxes Fa and Fb of the two species may differ. 

The migration of a surface atom is modeled as a 
nearest-neighbor hopping process with the rate Rd = 
ko cxp(— i^c/fc^T), where fcg = 10"'^^ Hz is an adatom vi- 
bration frequency, Ed is the hopping barrier, T is the 
substrate temperature and ks is Boltzmann's constant. 
The hopping barrier is the sum of a term from the sub- 
strate Esuh and a contribution from each lateral nearest 
neighbor E^. Both contributions depend on local compo- 
sition: For each term we have the four possibilities AA, 
AB, BA and BB. The hopping barrier of an atom X (of 
type A or B) is then 

E^= Y: (rr^E^ul+nrEn^ (3) 

Y=A,B 

where is the hopping barrier for a free X adatom 

on a substrate atom Y, Uq is equal to one if a substrate 
atom is of type Y and zero otherwise, n^^ is the num- 
ber of nearest-neighbor X-Y pairs, and E^^ is the cor- 
responding contribution to the barrier (symmetric in X 
and Y). Lateral interactions between impurity atoms are 
neglected (£;,f^ = 0). 

In the simulations reported in this paper we used 
Ef^i = 0.8 eV, EfX = 0.1 eV, Eg^ = 1.0 eV, Eg^ = 
0.1 eV, and the substrate temperature T = 500 K. The 
low values of E^f^ and E^^ ensure that atoms deposited 
on top of an impurity instantaneously descend to the 
substrate. Growth and impurity fluxes were varied in 
the interval ranging from 0.00025 ML/s to 0.25 ML/s. 
The system sizes ranged from 300x300 to 500x500. 
The nearest-neighbor coupling E^^ between A-atoms 
controls^^ the size of the critical nucleus i*. It may be 
stronger {E^^ > E^^) or weaker {E^^ < E^^) than 
the coupling to the impurities. In equilibrium at low 
temperatures, the former case leads to the formation of 
islands composed inside mainly of A atoms with B atoms 
bounded near the edges, while in the latter case it is en- 
ergetically more favorable when B atoms are inside the 
islands. 

However, our simulations show that growth leads to in- 
termixing of A and B atoms in both cases, -E^"^ > E:^^ 
and E^"^ < E^^ . Thus the energetic bias favoring seg- 
regation is not sufficient to obtain configurations with 
impurities mostly at island edges (decorated islands) . To 
achieve this, we have to introduce an additional thermally 
activated process, which allows an A atom approaching 
an island to exchange with an impurity covering the is- 
land edge. A similar process was introduced previously*^ 
in the context of homoepitaxy on Si(OOl) with prede- 



2 



posited hydrogen. In that work, an A-atom was allowed 
to exchange with an impurity provided the A-atom was 
not bonded to another A-atom at a nearest-neighbor site. 
In our case this modification turned out not to be suSi- 
cient, since impurities were still found to be progressively 
trapped inside islands during growth. We therefore allow 
the exchange of an A atom with an impurity also when 
it has a single bond to another A atom in a nearest- 
neighbor position. Using this rule, which is analogous 
to the exchange process invoked in the case of three- 
dimensional growth with surfactants,*'^® we obtain well 
decorated islands with impurities floating on the island 
edges during growth (Fig. lb, and c), see Section III. 

In principle, the rate of the exchange process could 
depend on the number of nearest-neighbor bonds (zero 
or one) of the A-atom. We observed that the differ- 
ence of both rates is not crucial provided both processes 
are active. The rates of these processes are taken as 
^ex = ko exp{—Eex/kBT), where Eg^ are the correspond- 
ing activation barriers. For both processes there is a 
maximum activation barrier above which the decorated 
geometry is not observed. In the following, the barriers 
for both types of exchange are for simplicity set to be 
equal. 



small E^^ = 0.1 cV, it is only partial for E^^ = 0.2 eV, 
and it becomes perfect for E^^ = 0.4 eV (cf. Figure 1). 
Hence, in order to get decorated islands, the barrier E^^ 
has to be larger than a minimal value. A simple detailed 
balance argument shows that the fraction /o of uncov- 
ered edge sites is given by 

/o = (l + 0Be^""/'=-^)-\ (4) 

and thus at T = 500 K a barrier of E^^ > 0.2 eV is 
required. As we shall see in Section V, the condition /o ^ 
1 also implies that the diffusion of j4-atoms is slowed 
down considerably by the impurities. 

Another important parameter determining decoration 
of island edges is the exchange barrier Ee^ that was in 
our simulations varied from 0.8 eV to 2 eV. For a small 
value of -Ecx, impurities are driven toward island edges, 
whereas for large i^ex, the exchange process is not active 
and impurities are often incorporated inside the islands. 
The impurities were observed to be floating on the island 
edges for smaU E^^, both for E^"^ > E^^ (Fig. lb) and 
for E^"^ < E^^ (Fig. Ic). Fig. Id shows a configuration 
for Sox = 2 eV in the case E^^ > E^^ . The surface 
morphology is similar to the one observed for E^^ < 
E^^ (not shown), where quite a regular checkerboard 



III. ISLAND MORPHOLOGY 

In Fig. 1 we show examples of typical configurations 
with the same partial coverage of both species 9a — 6b = 
0.1 ML (i.e., the total coverage 6 = 9a + 9b = 0.2 ML) 
obtained by codeposition of adatoms and impurities with 
fluxes Fa = Fb = 0.004 ML/s. Figures la, lb, Ic il- 
lustrate the effect of varying the relation between E^"^ 
and E^^ for E^^ = 0.3 eV. Several features can be 
identified. First, the island density increases with E^^ . 
This can be explained by the observation (cf. Fig. Ic 
for E^^ = 0.4 eV) that for larger E^^ , free A-atoms 
start to be captured by impurities and many small is- 
lands containing impurities and a few A atoms appear 
on the surface in addition to already existing decorated 
islands. These small islands act as nuclcation centers that 
lead to the increase of the island density. For large E^^ , 
almost all impurities will capture an A adatom for a cer- 
tain time. As we argue in Section V, this effect causes 
reduction of adatom mobility by the impurities. 

Whereas the island density increases, the density of 
free impurities decreases with increasing E^^ due to (i) 
the stronger A-B bond favoring the binding of impurities 
at island edges and (ii) the increase of the island density 
that leads to smaller islands with more perimeter sites. 
(Notice that we do not obtain decorated islands for simul- 
taneous deposition at very large fluxes or at very early 
stages of growth since there are not enough impurities 
available to cover all perimeter sites.) 

The degree of edge decoration also strongly depends on 
the value of E^^ . Edge decoration is not observed for 




(c) (d) 



FIG. 1. Examples of configurations for the total coverage 
= 0.2 ML obtained by simultaneous deposition with flux 
Fa = Fb = 0.004 ML/s, E^^ = 0.3 eV and different energy 
barriers E^'^ and Ec^: (a) E^'^ = 0.1 eV, Sox = 1 eV, (b) 
= 0.2 eV, Sex = 1 eV, (c) = 0.4 eV, E^^ = 1 eV, 
(d) ^ = 0.2 eV, Sex = 2 eV. We show only 50 x 50 sections 
of larger simulation boxes. 
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structure is produced with almost no free impurities on 
the surface. Thus in order to obtain decorated islands, 
i^ex has to be lower than a threshold value which in the 
present case is about 1.2 eV. 

A further remarkable feature of the configurations dis- 
played in Figure 1(b) and (c) is that the compact square 
island shape is maintained as the island density increases. 
In fact, careful inspection shows that the kink density on 
the well decorated edges is smaller than when the decora- 
tion is incomplete. This reflects the enhancement of edge 
diffusion by the impurities: The energy barrier for an 
adatom moving along the decorated step edge within the 
impurity layer is E^x which, under the conditions of Fig- 
ure 1(b), (c), is smaller than both the barrier + E^^ 
for diffusion along an uncovered step and the barrier 
-^sub + edge diffusion on the outside of the im- 

purity layer. Clearly this is true only if the exchange of 
singly bonded ^-atoms is allowed, which underlines the 
importance of this type of the exchange process. 

IV. ISLAND DENSITY SCALING 
A. Simultaneous growth 

In the previous section, we described rjualitativcly the 
growth morphology for one fixed value of the deposition 
flux. Here we present results for the behavior of the island 
density as a function of flux F and coverage 6, and 
discuss its dependence on the kinetic parameters. Results 
for each set of parameters were obtained by averaging 
over several independent simulation runs. 

In the presence of impurities, islands are composed of 
both A and B atoms. We define the size of an island 
as the number of A-atoms in a connected cluster of A- 
atoms forming the island. This definition is appropriate 
for growth with impurities segregating on the edges of 
the islands. However, visual inspection of configurations 
showed that for E'ex > 1-2, there exist islands contain- 
ing several mutually disconnected clusters of v4-atoms. 
Hence, our definition cannot be applied straightforwardly 
for large Ef.^. In the following, we restrict ourselves to 
situations where the intermixing inside the islands is neg- 
ligible. Simulations for larger values of i?cx indicate that 
the flux dependence of the island density flattens (the 
exponent x in (1) decreases), but due to the ambigu- 
ity in the definition of the island density in presence of 
intermixing, we did not attempt to assess the physical 
significance of this observation. 

1. Flux dependence 

Fig. 2 shows the island density iV as a function of the 
adatom fiux Fa for several coverages 9 and different en- 
ergy barriers E^^ and E^y^ (inset). The energy barrier 
E^-^ = 0.3 eV is fixed, and the impurities and adatoms 



are codeposited with the same flux Fb = Fa- For 
comparison we also show data for homoepitaxial growth 
without impurities at two coverages 6 = 0.05 ML and 
e = 0.1 ML. We can see that for E^^ = 0.2 eV and 
i?ex = 1 eV, the island density is quite close to the corre- 
sponding value in homoepitaxy. With increasing interac- 
tion energy between adatoms and impurities, the island 
density dramatically increases, but the exponent x in the 
power law relation (1) between flux and the island den- 
sity remains nearly unchanged. For example, we find 
X « 0.54 for E^^ = 0.2, x « 0.45 for E^^ = 0.4, and 
X ~ 0.54 for homoepitaxial growth, which means that the 
effective critical nucleus size is i* w 2 in this range of pa- 
rameters. According to Kandel's rate equation theory,^^ 
the scaling exponent should then become x ~ 0.8 in the 
presence of strong barriers to attachment. In our model 
this is not observed, because the bonding of the adatoms 
to the impurity-covered edges keeps them near the edge 
long enough for an exchange to occur. ^^'^'^ 

The inset of Fig. 2 shows that the island density is 
further increased if the exchange barrier £^ex is set to a 
larger value i^ex = 1-2 eV. The data for E^x = 1-2 eV and 
E^^ = 0.4 indicate a slight decrease of the exponent x- 

Figure 3 shows results obtained by varying the ratio 
Fb / Fa of impurity to adatom flux. In one set of simu- 
lations, using E^^ = 0.2 eV, the impurity flux was kept 
constant at Fb = 0.016 ML/s, while the adatom flux Fa 
was varied. For large fluxes Fa > Fb, the island density 



FA=Fa. E„ =0.3 eV 
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FIG. 2. Averaged island density as a function of flux Fa for 
several values of the total coverage 9 = 6a + Ob and different 
energy barriers: E^^ = 0.2 cV - open symbols, E^^ — 0.4 eV 
- filled symbols. The adatom interaction energy Et^ = 0.3 eV 
and the exchange barrier _Ecx = 1 eV are fixed, and the impu- 
rity flux Fb = Fa- The behavior in the absence of impurities 
(homoepitaxy, Fb = 0) is shown for comparison. Inset: The 
effect of the exchange barrier Box on the island density at 
coverage d = 0.2 ML. Circles represent the data from the 
main flgure for Bex = 1 eV, squares the corresponding data 
for Sex = 1.2 eV. 
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E„ =0.3 eV, E,^= 1.0 eV, codep. 9^=0.05 ML 



F.=0.016ML/s 




^ ^ V-""" 
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Fa 

FIG. 3. Dependence of the island density on flux Fa at 
constant Fb = 0.016 ML/s (open squares) compared with 
the situation Fb = Fa (open diamonds), and with homoepi- 
taxy (pluses), E^^ = 0.2 eV. Inset: The flux dependence at 
two constant ratios of fluxes, Fb/Fa ~ 2 (filled squares) and 
Fb/Fa = 1 (filled triangles), E^'^ = 0.4 eV compared with 
homoepitaxy (pluses). 

is seen to approach the data obtained for homoepitaxy, 
indicating that the impurities have no effect, while for 
small fluxes Fa < Fb, the flux dependence is described 
by an effective po-wer la"w N ~ "with x' ~ 0-36 < x- 
An interpretation of this behavior "will be given at the 
end of Section VB2 

In the second set of simulations shown in Figure 3, 
which were carried out using E^^ — 0.4 cV, both fluxes 
were varied keeping the ratio Fb /Fa = 2 constant. This 
is seen to further increase the island density without 
changing the flux dependence. In this sense, an increase 
in the coverage of impurities (by increasing Fb) is equiv- 
alent to increasing their effectiveness through an increase 
of the bond energy i?^^. A quantitative formulation of 
this statement will be given in Section V. 

2. Coverage dependence 

A new feature in comparison with homoepitaxy is a 
stronger coverage dependence of the island density. This 
is seen in Fig. 2 for both weak {E^^ = 0.2 eV) and 
strong {E^^ = 0.4 eV) interaction with impurities, but 
it is more pronounced for strong interaction, in particular 
at larger fluxes. We followed the coverage dependence in 
more detail for fixed flux, and in addition to the island 
density we also measured the density of free adatoms n. 
The results obtained at a medium flux Fa = 0.004 ML/s 
are compared with homoepitaxy in Fig. 4. Both in im- 
pure growth and in homoepitaxy the island density shows 
an initial regime of rapid increase followed by a "satura- 
tion" regime in which it increases much more slowly with 
coverage. However, the residual coverage dependence in 



F,=Fg=0.004 ML7S, E„ =0.3 eV 
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Coverage 9^ 

FIG. 4. Density of islands (solid lines) and free adatoms 
(dashed lines) as a function of coverage 6a for homoepi- 
taxy (no symbols) and for two different energy barriers 
E^^ = 0.2 eV (diamonds) and = 0.4 eV (circles). 

Adatoms and impurities are codeposited at the same flux, 
Fa = Fb= 0.004 ML/s. 

the saturation regime is stronger in the presence of impu- 
rities, and furthermore the onset of the saturation regime 
is delayed as the interaction between adatoms and impu- 
rities increases. A quantitative description of this effect 
will be provided in Section VB2. 

The density of adatoms exhibits a completely differ- 
ent behavior as compared to homoepitaxy. We observe 
that for gro"wth with impurities, the adatom density is 
comparable with the island density up to the coverage 
0A = 0.1 ML, whereas in homoepitaxy the adatom den- 
sity rapidly decreases after reaching a maximum at the 
beginning of the saturation regime (cf. Section VA). 
Other surprising features are the power-law increase n ~ 
Oa"^^ observed over almost two decades in the case of 
strongly interacting impurities, and the weak oscillations 
of the adatom density for coverages 9a > 0.01 (cf. Fig. 4). 
We will return to the behavior of the adatom density in 
Section VB2. 



3. Next-nearest-neighbor interaction 

A modification of the model in which the barrier for 
diffusion E^ contains an additional contribution from 
each lateral next-nearest neighbor of the opposite type 
was also studied. This implies that a term n2^E^^ is 
added to the right hand side of Eq. (3). Here ^ is 
the number of next-nearest neighbors of the type oppo- 
site to the atom under consideration, and E^^f is the 
corresponding contribution to the activation barrier. We 
do not consider ncxt-nearest-neighbor contributions from 
pairs of particles of the same type. For simplicity, the 
new parameter E^^ is set equal to E^'^. 
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FIG. 5. Examples of configurations for tlic total coverage 
= 0.2 ML obtained by simultaneous deposition with flux 
Fa = Fb = 0.004 ML/s, E^^ = 0.3 eV, Sex = 1 eV in a 
modified model with next-nearest-neighbor interaction: (a) 
E^^ = E^^^ = 0.1 eV, (b) = E^^ = 0.4 eV. We show 
only 50 x 50 sections of larger simulation boxes. 

Our motivation for introducing the additional term is 
a desire to stiidy an improvement in the decoration of 
island edges by impurities, and the resulting decrease of 
the fraction /o of uncovered edge sites. The additional 
interaction also enhances the probability of nucleation 
around impurities because the number of sites at which 
an adatom can be captured is considerably higher. The 
configuration shown in Fig. 5a demonstrates that now we 
obtain almost perfect decoration also for E^^ = E^^ = 
0.1 eV. Fig. 5b illustrates the nucleation of small islands. 

For E^^ = E^^ =0.1 eV, the island density is nearly 
the same as for homoepitaxy for all fluxes studied, and 
decoration is perfect provided there is sufficient amount 
of impurities available. This shows that also in the pres- 
ence of next-nearest-neighbor interactions the decorated 




FIG. 6. Examples of configurations deposited with fiux 
Fa = 0.004 ML/s at the coverage 61^ = 0.1 ML after pre- 
deposition of Ob = 0.1, using parameters: E^"^ = 0.3 eV, 
Sex = 1 eV, (a) = E^^ = 0.2 eV, (b) 

= E^J^ = 0.4 eV. We show only 50 x 50 sections of 
larger simulation boxes. 



edges are unable to block efficiently the attachment of 
adatoms. For larger values of E^^ = E^^ , the addi- 
tional interaction causes an increase of the island density 
and a decrease of the scaling exponent x- For example, 
the effective value for E^^ = E^^ = 0.2 eV is x « 0.42 
and for E^^ = E^^ = 0.4 eV it drops to x « 0.3. This 
suggests that the nucleation of small islands described 
above eff'ectively lowers the size i* of the critical nucleus. 

B. Predeposition of impurities 

We performed simulations with predeposition of impu- 
rities for the same set of parameters as for codeposition. 
In order to obtain a morphology with island edges deco- 
rated by impurities, we need an appropriate value of i?ex- 
Complete decoration also requires a sufficient amount of 
impurities available on tlw^ surface. The data presented 
are for a predeposited coverage Ob = 0.1 ML. Examples 
of morphologies for E^^ = 0.2 eV and E^^ = 0.4 eV are 
shown in Fig. 6 and look qualitatively similar to Figures 1 
b and c. 

The i^yi-dependence of the island density is compared 
with the results for codeposition in Fig. 7. The island 
densities in the predeposition regime are slightly higher 
than for codeposition. This is qualitatively plausible, 
since the predeposited impurities are present on the sur- 
face throughout the deposition process and hence their 
effect on growth accumulates over time. The correspond- 
ing curves are shifted by a factor independent of the flux. 
The slope remains the same as for codeposition. The dif- 
ference from codeposition is that there is no appreciable 



E„ =0.3 eV, E„= 1 .0 eV 
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FIG. 7. Comparison of averaged island density as a func- 
tion of flux Fa for codeposition (circles) and predeposition 
(triangles) for different energy barriers: E^'^ = 0.2 eV - open 
symbols, E:^^ = 0.4 eV - filled symbols. The predeposited 
coverage is 9b = 0.1 ML. The adatom interaction energy 
E^^ — 0.3 eV and the exchange barrier Sox = 1 eV are fixed. 
The behavior for homoepitaxy is shown for comparison. 
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F^=0.004 MUs, 9j=0.1 , E„ =0.3 eV 
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FIG. 8. Density of islands (solid lines) and free adatoms 

(dashed lines) as a function of coverage Qa for hornocpitaxy 
(no symbols) and for prcdeposition with two different energy 



barriers = 0.2 eV (diamonds) and 



0.4 eV (cir- 



cles). Flux of adatoms is Fa = 0.004 ML/s, the predeposited 
coverage of impurities is 9b = 0.1 ML. 

coverage dependence for Oa larger than 0.05 ML, except 
at larger fluxes in the case of strongly interacting impu- 
rities. 

The detailed coverage dependence for a fixed flux is 
shown in Fig. 8. We see that the behavior of both the is- 
land density and the adatom density is qualitatively sim- 
ilar to homoepitaxy. The island density saturates and at 
the same time the adatom density starts to decrease. The 
impurities only cause a shift of the crossover to the satu- 
ration regime to a higher coverage, the shift being larger 
for stronger interaction between adatoms and impurities. 
We shall return to this eS^ect in Section VB 1. 
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(6) 



has been introduced. 

The first conclusion that can be drawn from Eq. (5) is 
that predeposited impurities, 6b = const., significantly 
affect the adatom diffusion only if 9b ^ !/(/'. In the 
case of codeposition 0b = Fst and D becomes time- or 
coverage-dependent. It is then useful to rewrite (5) in 
terms of the coverage 6a of ^-atoms as 



D 



D 



1 + 6aI6* ' 



with the characteristic coverage 

6* = 



(7) 



(8) 



For coverages 9a > 9* the impurities begin to signifi- 
cantly aH'ect the adatom mobility. The expression (8) 
quantifies the statement made above in Section IV A 1 
that an increase of the flux ratio Fb /Fa in codeposition 
is eqiiivalent to an increase of the A-B binding energy 
E^^. For Fa = Fb and T = 500 K, we have 6* « 0.01 
for E^^ = 0.2 eV and 6* « 10"* for E^'^ = 0.4 eV. In 
the following, these two sets of parameters will be referred 
to as the case of weak and strong impurities, respectively. 



A. Pure growth 

We proceed by combining (5) with the simplest ana- 
lytic model of nucleation, consisting of two coupled rate 
equations for the island density and the adatom den- 
sity n. In the absence of impurities, the equations for a 
critical island size i* — 1 read^^'^^'^'^ 



djTh 

— =Fa- 4£>n(2n -|- N) 



(9) 



V. RATE EQUATION THEORY 

In this section, we develop a simple rate equation ap- 
proach to explain, at least qualitatively, the main impu- 
rity effects on the island density which were presented in 
Section IV. Our basic assumption is that the impurities 
affect the growth process only by slowing down the diffu- 
sion of adatom,s. To obtain a simple analytic expression 
for the effective adatom diffusion coefficient D(9b) in the 
presence of an impurity coverage 9b, we further replace 
the (mobile) impurities by static traps with binding en- 
ergy E^^ . Then standard results for diffusion in random 
media yield^^ 



D{9b) 



D 



D 



(5) 



where D = fcoe^^i^>b/'^B^ jg diffusion coefficient of a 
single adatom on the clean substrate, and the abbrevia- 
tion 



— = ADn^. 
dt 



(10) 



The main features of the solution of (9,10) with initial 
condition n = = can be described as follows (sec 
the article by Tang^° for a lucid presentation): In the 
early time regime the adatom density increases linearly 
by deposition, n « F^t = 6a-, and accordingly the island 
density grows as 



N ^(A/?,){D/Fa)9\ 



(11) 



In the late time regime the adatoms are mainly captured 
by preexisting islands. This implies that n w Fa/DN <C 
N and the island density grows more slowly, as 



N « {FA/l2DfN^f, 



(12) 



"1/3 

while the adatom density decreases as n ~ 6*^ . The 
transition between the two regimes occurs at a coverage 
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^1 ~ {FA/Dfl\ 



(13) 



Keeping the eovcragc fixed while increasing the flux 
therefore takes the system from the late time regime, 
where N ~ F^/^, into the early time regime with N <^ 
F~^, with a maximum in the island density attained at 
a critical flux F" ~ 00"^. 

To generalize these estimates to the case i* > 1, we 
replace the nucleation equation (10) by^^'^* 



Dn 



(14) 



Then the early time behavior becomes N ~ {D / F)0\~^^ , 
while in the late time regime 



+2) 



(15) 



in agreement with the expression (2) for the scaling ex- 
ponent X- The transition coverage 9i is estimated by 
matching the two behaviors, which yields 



~ (F^/£»)2/(i*+3). 



(16) 



For i* = 1 the coverage dependence of the densities 
of islands and adatoms observed in microscopic simula- 
tions is in accordance with the rate equation theory.^°'^^ 
In the reversible case i* > 1, the simple rate equations 
are quantitatively inappropriate, though the key qualita- 
tive features - the existence of an early time regime of a 
rapid increase of the island density, followed by a "pre- 
coalescence saturation regime" with little change in iV - 



remain 



22 



B. Impure growth 

1. Predeposition 

The effect of predeposited impurities is obtained sim- 
ply by replacing D by the constant expression (5) for D 
in (15) and (16). Consequently the island density in the 
late time regime increases by a factor {1 + 9B(t>y ^'■^ 
which is independent of flux or ^-coverage, and the on- 
set of saturation is delayed by a factor (1 -|-^b<^)^/('*+3\ 
This is in qualitative agreement with the coverage de- 
pendence of densities of islands and adatoms displayed 
in Figure 8, which shows the same overall behavior as in 
the homoepitaxial case, only shifted to larger coverages 
and higher densities. Quantitatively, the numerically ob- 
served increase in the island density is consistent with the 
factor (1 + ^B^)* if the size of the critical nucleus 

is set to i* = 1. On the other hand, if the critical nucleus 
size is assumed to be z* = 2 as suggested by the numer- 
ical value of Xi then the theory is seen to overestimate 
the increase in the island density and underestimate the 
delay of the onset of saturation. 

As was mentioned in Section V A, the island density at 
fixed coverage 6a shows a maximum at a critical flux F'^, 



which is determined by setting the saturation coverage 
equal to 9a- For predeposited impurities with 9b4> ^ ^ 
this is given by 



Fl ^ {D/c^9b)9'{+'^/'. 



(17) 



Further discussion of predeposition in relation to code- 
position will be provided below. 



2. Codeposition 

In the case of codeposition the situation is richer due 
to the coverage dependence of D. First, the cases 9i < 9* 
and 6*1 > 9* have to be distinguished. In the first case the 
impurities only affect the late time regime. For 9a ^ 9* 
we can approximate (7) hy D D9* /9a- Inserting this 
into the nucleation equation (14) and setting n « F/DN 
we obtain the expression 



TV w {Fa/D9*) 



i'/ii'+2) 



{Oa) 



(i- + l)/(i*+2) 



(18) 



which replaces (15) in the pure case. It can be seen that 
the scaling of N with flirx Fa remains the same, i.e. the 
exponent x is not affected. The impurities increase the 
island density by a factor (0^/0*)'*/^'*+^) which, in con- 
trast to the case of predeposition, is coverage dependent. 
This is qualitatively consistent with the residual coverage 
dependence of the island density seen in Figure 2, which 
is nearly absent in the corresponding predeposition data 
in Figure 7. In quantitative terms, however, the cover- 
age dependence in (18) is much stronger than what is 
observed in the simulations. The dependence of (18) on 
9* shows that increasing the ratio Fb / Fa at constant Fa 
and 6a will also increase the island density, in accordance 
with the simulation data shown in Figure 3. 

Noting that for codeposition 9b = {Fb/Fa)9a = 
9a/{(I>6*), Eq. (18) can be rewritten as, N 
{FA9B(t>/Dy'^^''+^^9]/'''^^\ which is identical to the ex- 
pression for predeposition with 9b4>'^ 1- A more careful 
calculation shows that the island density for predeposi- 
tion exceeds that for codeposition by a constant factor 
{i* + l)i/("+2), if systems of the same impurity coverage 
6b are compared. This is qualitatively consistent with 
Figure 7, though the simulation data indicate that the 
factor is larger for weak impurities than for strong ones. 

A surprising consequence of the rate equations with 
coverage-dependent diffusion is that the adatom density 
increases with coverage in the late time regime. This 
follows from balancing the deposition term on the right 
hand side of (9) against the island capture term ADNn « 
AD{6*/9A)Nn, which yields 



(FA/i?r)2/(«*+2)^y(^*+2). 



(19) 



In fact, when 6i < 9* the adatom density is a nonmono- 
tonic function of coverage: It increases as n ^ for 
9a<9i, decreases as n ~ 6^^^^''^^^ for 6i<6a<9*, 
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F.=0.004 MUs 



6=0.10 ML 




FIG. 9. Time evolution of the adatom (open symbols 
and the dashed line) and island (full symbols and the 
solid line) densities obtained from a numerical solution 
of the rate equations for i* — 1 [Eqs. (9) and (10)] 
with the coverage-dependent diffusion constant (5). The 
parameters were chosen for comparison with the simula- 
tion data in Figure 4; Flux Fa = Fb = 0.004 ML/s, 
D = (A;o/4)exp(-S,4b/fesT) « 2.15 x 10*8'^, and 
E^'^ = 0.2 eV (diamonds) and E^'^ = 0.4 eV (circles). The 
agreement between the homoepitaxial island density (solid 
line) and the adatom density for weak impurities (open cir- 
cles) is coincidental. 

and increases again according to (19) for 6a > 6* . This is 
illustrated in Figure 9, which is reminiscent of the simula- 
tion data for the adatom density in Figure 4. However, in 
contrast to the simulations, the late time adatom density 
given by (19) is small compared to the island density (18), 
and neither the intermediate scaling regime n ^ {dAf'"^^ 
nor the oscillations of n seen in Figure 4 are reproduced 
by the rate equations. Here and in the following figures 
we show results obtained by numerical integration of the 
rate equations (9,10) for i* = 1, with D replaced by D. 
This is sufficient for a qualitative comparison, and re- 
lieves us of the necessity to explicitly treat the growth 
dynamics of the intermediate unstable clusters. 

Consider next the case 9\ > 9*. For coverages in 
the early time regime which satisfy 6a ^ 6*, we set 
D K, DO* /Oa and obtain, using n k, 9a, the early time 
behavior N ~ {D6* / Fa)6\+^ . The onset of the satura- 
tion regime then occurs at a coverage 

9i ~ {FA/De* f/^^+'''> (20) 

which exceeds the corresponding expression (16) for the 
pure system by a factor (^i/^*)^/^'*"'"^^ Thus both pre- 
deposited and codeposited impurities delay the onset of 
the saturation regime. 

The critical flux at which the island density attains a 
maximum is now given by 



FIG. 10. Rate equation predictions for the flux dependence 

of the island density in the case of codeposition with Fa = Fb 
at the total coverage 6 = 0.1 ML/s. The parameters were 
chosen as in Figure 9. For comparison with Figure 6, data 
for predeposition with 9b = 9a = 0.05 ML/s in the case of 
strong impurities are shown as well. 

F% ~ (D/0)0f +1)/^ (21) 

which is seen to become identical to the predeposition 
expression (17) by setting 9b = {Fb/Fa)9a- A quanti- 
tative evaluation using the parameters of the simulations 



codeposition, 


8,=0.05 lUIL 


Fb=0.016ML/S 




— - F,=F3 




homoepitaxy 








^ ' " ^ 


10^ 10'^ io' 


io' 












/ lO"' 









10* 10"' 10"^ 10"' 



Fa 

FIG. 11. This figure illustrates the rate equation 

prediction for codeposition with constant impurity flux 
Fb ~ 0.016 ML/s (solid line) and should be compared to 
Figure 3. The predictions for the case Fa = Fb (dashed 
line) and for homoexpitaxy (dotted line) are also shown. 
E^^ = 0.2 eV, other parameters were chosen as in Figure 9. 
Inset: The rate equation prediction for codeposition with con- 
stant flux ratio Fb/Fa = 2 (long-dashed line) compared to 
the case Fa = Fb (dashed line) and to homoexpitaxy (dotted 
line), E^'^ = 0.4 eV. 
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shows that the critical flux is beyond the range of sim- 
ulated fluxes for the case of weak impurities, while it 
should be observable for strong impurities. This is il- 
lustrated in Figure 10. It therefore appears natural to 
identify the saturation of the flux dependence of the is- 
land density found in the simulations (see Figures 2 and 
7) with the maximum predicted by the rate equations, 
although it should be emphasized that the simulations 
show no clear evidence of a decrease of N beyond the 
plateau. 

Finally, we address the effect of changing the adatom 
flux Fa while keeping the impurity flux Fb constant. 
At fixed 6a this implies a crossover from 9a <C 9* 
for Fa > Fb<I)9a to 9a > 9* for Fa < -Fs^fA- 
In the high flux regime the island density A'' is unaf- 
fected by the impurities, while in the low flux regime 
N - {Fa/DY{9a/9*Y - {Fb<P/DY becomes indepen- 
dent of Fa, since 9* ~ Fa- This is illustrated in Fig- 
ure 11, which should be compared to Figure 3. The 
behavior is qualitatively similar, however instead of a 
plateau at low fluxes the MC simulations show a sec- 
ond scaling regime N ^ with a nonzero exponent 
X' < X- 

3. The effect of edge decoration 

In the formulation of the rate equations we have as- 
sumed that all deposited impurities contribute to the 
impurity coverage in the expression (5) for the effective 
diffusion coefficient, thus neglecting the fact that a cer- 
tain fraction of impurities is bound at the island edges. 
This assumption is self-consistent only if the density rig 
of edge sites, as predicted by the rate equations, is small 
compared to the deposited coverage 9b oi impurities at 
all times. 

For compact islands the density of edge sites is of the 
order of N\/~A, where A ~ {9a — n)/N is the area of an 
island, and hence 

ne ~ ^/N{9a - n). (22) 

In the saturation regime n, N <^ 9a, thus Ue ~ \/9aN. 
In the early time regime n w to leading order, and the 
density of edge sites is determined by tlic^ iicxt-to-lcading 
correction. Analysis of the rate equations shows that 
He ~ AT both for 9a < 0* and for 6a > 9*, which imphes 
that the island size does not increase with coverage in 
this regime. 

Using the estimates for rig, the importance of the im- 
purities bound at the edges can be worked out for spe- 
cific cases. Since evidently Ug ^ 6a always, a sufficient 
condition for the irrelevance of edge decoration is that 
6b > 9a &t all times. This is true for predeposited im- 
purities up to a coverage 9a = 9b, and for codeposition 
with Fa < Fb- Among the situations treated earlier in 
this section, the only case where corrections due to edge 



decoration may be expected is deposition at fixed impu- 
rity flux Fb in the transition region Fb < Fa < Fb4>9a 
(see Figure 11). For fluxes Fa > Fb4>9a the impurities 
were seen to be irrelevant even if the full impurity cover- 
age contributes to slowing down the adatoms. Since the 
edge decoration decreases the impurity concentration on 
the terraces, it is likely that its only effect will be to shift 
the point where the island density becomes equal to its 
homoepitaxial value towards smaller deposition fluxes. 

VI. CONCLUSIONS 

A brief glance at the results presented in this paper 

may lead to disappointment: No dramatic change of the 
exponent x has been found contrary to expectations in 
any of the modiflcations of the simulation model. ^'^ How- 
ever, upon closer inspection our work reveals several non- 
trivial features that should not be overlooked in the large 
amount of numerical data. 

(i) We have established that a perfect decoration 
of island edges in our model requires a process of 
adatom exchange that is completely analogous to the 
exchange discussed for surfactant-mediated growth of 
semiconductors.^^ A feature worth remembering is the 
enhancement of diffusion of adatoms along island edges 
via the exchange process with impurities attached to 
these edges. This mechanism of smoothing on a one- 
dimensional substrate provides another perspective and 
possible interpretation of the smooth growth on a two- 
dimensional substrate in the; presence of a surfactant. 
Smoothing of island shapes should be experimentally ver- 
ifiable and may have practical implications. 

Given the fact that size and shape of islands (including 
the number of kinks at island edges) in the submonolayer 
regime of growth determine the developing surface mor- 
phology (cf. most of the experimental papers cited below, 
and in particular Ref. 10), a possibility to control both of 
them by adding impurities seems very attractive. 

(ii) From a theoretical perspective, we have been able 
to obtain insight into (and even semi-quantitative agree- 
ment with) the simulation results, using rather simple 
rate equation theory. The surprising resistance to change 
of the exponent % can be understood within this theory, 
as well as other features of the simulations, such as the 
strikingly different behavior of the adatom density dur- 
ing impure growth with codeposition as compared to the 
case of homoepitaxy. 

Our research also leads to new questions. One would 
like to understand better the details of behavior observed, 
in particular the oscillations of the adatom density seen 
in simulations (Fig. 8). Another open question concerns 
"monovalent" impurities that can bond to island edges 
but do not bond to adatoms approaching these edges, and 
which therefore should bring about an efficient passiva- 
tion of the island boundaries. Preliminary simulations 
suggest that this effect alone does not bring about any 
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significant change in the value of the exponent %. 
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